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Deuteron photodisintegration with polarized photons at astrophysical energies 

G. Ramachandrar0 and S. P. Shilpashree 
Indian Institute of Astrophysics, Koramangala, Bangalore-560034, India 
(Dated: February 9, 2008) 

Following precise experimental studies at the Duke Free-Electron Laser Laboratory, we discuss 
photodisintegration of deuterons with 100% linearly polarized photons using a model independent 
theoretical approach taking together Ml and El amplitudes simultaneously. The isoscalar Mis 
contribution is also taken exactly into account. From the existing experimental measurement on 
doubly polarized thermal neutron capture, it is seen that the isoscalar Mlg contribution could be of 
the same order of magnitude as the experimentally measured cross sections at energies relevant to Big 
Bang Nucleosynthesis (BBN). Therefore appropriate measurements on deuteron photodisintegration 
are suggested to empirically determine the Alls contribution at astrophysical energies. 

PACS numbers: 23.20.Gq, 29.25.Pj,26.35.+c, 24.70.+S, 21.45.+V, 25.20.-x 



The photodisintegration of the deuteron and its in- 
verse reaction viz., n-p fusion in the neutron energy range 
of order 10 to 10^ keV is of considerable interest to as- 
trophysics. It is important to the nucleosynthesis sce- 
narios fll fr om the Big Bang to stellar evolution under 
various conditions. The earliest estimates of the reac- 
tion rates by Fowler, Caughlan and Zimmerman (FCZ) 
used theoretical calculations Q of deuteron photo- 
disintegration normalized to the then available thermal 
neutron radiative capture cross section measurements 
P|. In a comprehensive evaluation of the reaction rates 
and uncertainties in 1993, Smith, Kawano and Malaney 
have pointed out : ^^With a binding energy of 2.22 
MeV, deuterium is the most fragile of the primordial 
isotopes: it is rapidly destroyed in stellar interiors.... 
Given that significant quantities of deuterium can only 
be produced during primordial nucleosynthesis, detection 
of deuterium provides important evidence in favor of the 
big bang model... Given the range of D/H observed in 
the interstellar medium, it is difficult to directly deter- 
mine a lower limit. . . a determination of the upper limit 
is plagued by uncertainties arising from chemical evolu- 
tion effects. The ratio of the primordial abundance of 
deuterium to that observed today could be any where be- 
tween 1 and 50^' . 

Laboratory measurements and decisive developments 
in astronomical observations go hand in hand to remove 
crucial ambiguities in nuclear physics input parameters 
and sharpen theoretical predictions in the astrophysical 
context. Although laboratory measurements with ther- 
mal neutrons date back to 1936 by Fermi and collabora- 
tors 1^ it has not been possible for a long time to measure 
the cross section at astrophysical energies due to the ten- 
dency of the neutrons to thermalize at low energies. The 
first cross section measurements between 20 keV and 64 
keV have been reported in 1995 by Suzuki et ai, and 
subsequently, by Nagai et al., at 550 keV. Buries and 
Tytler |^ measured the deuterium abundance in high- 
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red-shift-hydrogen clouds, where (it may be expected 
that) almost none of the deuterium could have been de- 
stroyed subsequent to the primordial stage. However, in 
a re-examination of the estimates of the uncertainties in 
1999 to sharpen the predictions of big bang nucleosyn- 
thesis (BBN), Buries, NoUett, Truran and Turner 
have observed : " Our method breaks down for the pro- 
cess p-\-n — > d-\-j. This is because of a near-complete lack 
of data at the energies relevant for BBN. The approach 
used for this reaction is a constrained theoretical model 
that is normalized to high precision thermal neutron cap- 
ture cross-section measurements^^ . The measured cross 
section of 334.2 ± 0.5m6 by Cox, Wynchank and Collie 
|lll | for thermal neutrons is considered as standard. 

The thermal neutron cross section has traditionally 
been interpretted in terms of a dominant isovector A/1 
amplitude for radiative capture from the initial ^Sq state 
of the n-p system in the continuum. Theoretical cal- 
culations [3 based on potential models led to a 10% 
discrepancy with the experimental measurements. Breit 
and Rustgi 13] proposed a polarized target-beam-test to 
detect the possibility of radiative capture from the initial 
^Sq state as well, which can take place through isoscalar 
A/1 and possibly also isoscalar E2 transitions. However, 
the surprising accuracy with which Riska and Brown [T^ 
explained the 10% descrepancy by including Meson Ex- 
change Current (MFC) contributions, set the trend for 
theoretical discussion in later years. It has been noted 
by Nagai et ai, ^ that the measured cross section is 
in a gree ment with the theoretical calculations by Sato et 
at, [13 including MFC's, isobar currents and pair cur- 
rents. They have also pointed out that "the theory is 
in good agreement with the cross section measured for 
neutrons above I4 MeV, but it deviates by about 15% 
from the measured cross section of the d{'y, n)p reaction 
by using the 7 ray of between 2. 5 and 2. 75 Me V ^3/, 
corresponding to neutron energies of 550 and 1080 keV^^ 
8J. Experimental studies on photodisintegration of the 
deuteron for photon energies from 2.62 MeV and above 
is well documented 17] . The cross section at 2.62 MeV is 
1.30±0.029m& which increases slowly to 2.430±0.17to5 at 
4.45 MeV and starts slowly decreasing with energy there- 
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after. The disintegration process is dominantly through 
El transitions leading to final triplet P-states of the n-p 
system in the continuum. Apart from the 15% discrep- 
ancy with the measured cross section noted by Nagai et 
al., 1^, the measured angular distribution and neutron 
polarization at photon energy of 2.75 MeV |0| and in the 
range 6 to 13 MeV ligl were found to be in disagreement 
with theoretical predictions which included the meson ex- 
change currents. Measurements of the analyzing power 
pof in p{n, j)d at neutron energies of 6.0 and 13.43 MeV 
were consistent with and theoretical calculations |2ll | 
showed that meson exchange currents produce a signifi- 
cant change but the effect is to move the theoretical curve 
to more negative values, thus making the discrepancy be- 
tween theory and experiment more pronounced. An ob- 
servable which is sensitive to the presence of isoscalar Ml 
and E2 transitions from the triplet 5-state is the circular 
polarization of the emitted radiation with initially polar- 
ized neutrons. The first measurement j22j to detect the 
presence of isoscalar amplitudes was not quite encourag- 
ing but a subsequent measurement |2^ yielded a value 
P-y = -(2.29 ± 0.9) X 10-3. An attempt ^ to explain 
the large measured value by introducing a six quark ad- 
mixture in the deuteron wave function led however to 
a disagreement with the well known deuteron magnetic 
moment. Later calculations "25l in the zero range ap- 
proximation and the wavefunction for a Reid soft core 
potential led to a theoretical prediction of the order 
of —1.1 X 10"'^ with an estimated accuracy of 25%. The 
measured value |2^ of P^ = — (1.5±0.3) x 10"'^ is in rea- 
sonable agreement with the theoretical calculation [25| . 
The importance of measuring the photon polarization 
with initially polarized neutrons incident on a polarized 
proton target has been pointed out |23| • When the initial 
preparation of the neutron and proton polarizations P{n) 
and P{p) are such that they are either opposite to each 
other or orthogonal to each other, the interference of the 
small isoscalar amplitudes with the large isovector am- 
plitude could substantially contribute to the observable 
photon polarization. 

Anticipating the experimental results of polarized ther- 
mal neutron capture by polarized protons by Miiller et 
al, the possibility of the initial ^S'l state contri- 

butions at thermal neutron energies was discussed using 
two different versions of effective field theory [2^ l30l | Al- 
though the measured value of (1.0 ± 2.5) x 10~^ for the 
7 anisotropy rj was not sufhciently sensitive to distin- 
guish between the two theoretical predictions, we may 
use equation (2) of Miiller et at, |23] to estimate the ra- 
tio R of the triplet to singlet capture cross sections to 
be 1.202 X 10~3. If we multiply R by the well-known 
cross section Tl], we get an estimate of 401.7 iib for the 
■^5*1 contribution to the cross section at thermal neutron 
energies. Quite surprisingly, this number is of the same 
order as the measured cross sections for capture at as- 
trophysical energies of 20, 40 and 64 keV [3. In fact, it 
is even larger by a factor of 10 than the measured cross 
section at 550 keV 0. This raises an open question as 



to what could possibly be the ratio R at astrophysical 
energies relevant to BBN. 

The influential paper of Buries, NoUett, Truran and 
Turner has inspired several theoretical [s^ 
as well as experimental [s^ studies. Since pho- 
todisintegration of the deuteron is well documented 
35] for photon energies of 2.62 MeV and above and 
is known to be dominated by El transitions leading 
to final triplet P-states in the n-p continuum, these 
studies were motivated towards the determination of 
the relative Ml and El contributions to the process 
at astrophysical energies. The experiment |33| was 
concerned with the measurement of the near threshold 
beam analyzing power using for the first time a laser 
based 7-ray source at 3.58 MeV. This was followed by 
measurements at seven 7-ray energies between 2.39 and 
4.05 MeV HI]. These measurements with 100% linearly 
polarized photons have been analyzed, making several 
simplifying assumptions viz., 

a) only I = 0,1 partial waves were considered in the final 
state due to the low energies involved, 

b) of the allowed two Ml and four El transitions, the 
isoscalar El leading to ^Pq is set to zero, 

c) the isoscalar Ml term leading to ^Si is neglected, 
using the traditional agruments for its supression, 

d) the three isovector El terms were combined to form 
a single P-wave amplitude, 

using the theoretical formalism [s^, where Ml and El 
contributions were calculated separately. 

The purpose of the present paper is to study (i(7, n)p 
theoretically, using a model independent formalism, 
without making any simplifying assumptions except that 
only the dipole transitions are considered with I = 0,1 
partial waves in the final state. Since the strength of the 
isovector Mly amplitude which is dominant at thermal 
neutron energies is known to decrease [isl l3ll Is^ by 
several orders of magnitude as energies relevant to 
BBN is approached and an estimate of 401.7 iib of the 
contribution of the isoscalar Mlg amplitude to the cross 
section at thermal neutron energies is seen to be of the 
same order of magnitude as the measured cross sections 
0,0 at energies relevant to BBN, it is not unreasonable 
to pay attention to the contribution of the isoscalar 
Mis amplitude at the energies of astrophysical interest. 
Moreover, spin observables are generally sensitive to 
the interference of a leading amplitude with other 
amplitudes which are not expected to be large. It is 
therefore appropriate to study the sensitivity of the 
beam analyzing powers to the isoscalar Ml^ amplitude 
leading to final ^Si state at astrophysical energies. 

We choose the linearly polarized photon momentum k 
in cm. frame to be along z-axis and the linear polariza- 
tion to be along x-axis of a right handed cartesian coordi- 
nate system and the neutron momentum p in cm. frame 
to have polar coordinates {p,9,4>), following [s^]- If the 
left and right circular states of photon polarization are 
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defined following Rose JST'I through = — /i^^, ^ = ±1, 
the above state of linear polarization may be represented 
by -^(u+i + u_i). We use natural units, h — c — 1. 
The unpolarized differential cross section for the reaction 
(i(7, n)p, in cm frame at energy E is given by 



dn 



1 E,,EpEd\p\ 



6 (27r£;)2 

I J2 rr[M(^)Mt(^)] 



^ Tr(T(A.)Tt(^)) 
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(1) 



where Tr denotes the trace or spur and T(^) denotes the 
on-energy-shell matrix for d{'y, n)p, when photons are in 
the polarized state u^. The cm. energies of the neutron, 
proton and deuteron are denoted respectively by E. 
and Ed- Following we express 



n ; Ep 



1 s+1 

M(m)=E E {S\s,l)■:F\s,^^)), 

s=0 A=|s-l| 



(2) 



in terms of irreducible tensor operators, S^{s, 1) of rank 
A in hadron spin space [39l | connecting the initial spin 
1 state of the deuteron with the final singlet and triplet 
states, s = 0, 1 of the n — p system in the continuum. 
Making use of the multipole expansion for u^e^''''" [s^l 
and expressing the continuum states of the n-p system in 
terms of partial waves, the irreducible tensor amplitudes, 
J-^{s, n) of rank A are given, in general, by 



completely characterized at any energy by the set of four 
irreducible tensor amplitudes T^{s, fj,), given by ©. But 
the contributing partial wave multipole amplitudes F^^J.j^ 
increase as the cm. energy increases. 

In the region of interest to BBN, we may restrict our- 
selves to only L = 1 and to / = 0, 1 partial waves as in 
|33l |. Then, we clearly have two All amplitudes viz., the 
isovector Ml^ leading to the final ^5*0 state, the isoscalar 
Mis leading to the final state and four El ampli- 
tudes viz., three isovector iJlj^^'-'^'^ leading to the final 
^Pj states and an isoscalar Elg leading to the final ^Pq 
state. In terms of these limited number of partial wave 
multipole amplitudes, the four irreducible tensor ampli- 
tudes ^l){s, /i) may explicitly be written as 

.^1(0,^) - -zA/l,/i(0,l,M)-A/3i?l,/i(l,l,/i), (6) 
= i£;i.(0)/0(l,l,Ai), (7) 



:EUl)fl{l, 1, fi) + iMlJliO, 1, /i)(8) 
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EU{2)f^{l,l,^^). 



(9) 



The £'1„(A) amplitudes contributing to the triplet irre- 
ducible tensor amplitudes JF^(l,/i) with A = 0,1,2 are 
related to the Ell aniplitudes with j = 0, 1, 2 through 



(10) 



EU{0) 




"13 5 " 






El^il) 




2 3-5 




EV=^ 

V 


El,{2) 




2 -3 1 







= ^EE E E W 

L = l i=0 3 = \l-s\ /=0,1 



l-\-s 



(3) 



where denote the orbital angular momentum and 
isospin in the final state, j denotes the conserved total 
angular momentum, L denotes the total angular mo- 
mentum of the photon and the shorthand notation \V\ 
stands for \/2L + 1. The partial wave multipole ampli- 
tudes fI^j^ depend only on cm. energy E, while the 

X C{l,L,\;mi,-^i,i^)Yi„i,{9,(l3), (4) 

take care of the angular dependence and also the depen- 
dence on photon polarization. The projection operators 
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^^-'} (5) 

assume either of the values 0,1 such that, if 7r+ = 1 im- 
plies TT^ = and vice versa. The F^J.j^ denotes electric 

2^-pole amplitudes, if 7r+ = 1 and magnetic 2^-pole am- 
plitudes, if TT" = 1. It may be noted that the reaction is 



The differential cross section relevant to [33, for 
d(7, n)p with linearly polarized photons is given, in cm. 
frame, by 



da 

dn 



TrMM\ 



where 



M = M(+l) + M(-1) 



(11) 



(12) 



Using known properties p3a | of the irreducible tensor op- 
erators and standard Racah algebra, we have 



da 27r2 r , 9 „ / 

— = [a + b snr 9(1 + cos 20) 

dil 6 



c cos ( 



(13) 



where 



[8|MU|2 + 24|M1«|2 + 36\Els\^ + 8|i;irT 
18|£;irT + 26|i;ir'l' - 16ReiEli=°Eli=^*) 
36i?e(£;ir'i?ir'*)]> (14) 



b = [9\Ell=^\^ + 21\Eli='^\'^ + 2'iRe{Eli=°Eli='^*) 
+ 5ARe{Eli-^Eli=^*) - 18\Els\^], (15) 
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and 

c = AV6Re[{2Eli=° + 3Eli=^ - 5Eli='^)Ml*] . (16) 

It is readily seen from H16(l that the third term c cos 9 in 
arises due to the interference of the Mlg amphtude 
with the Ely ampHtudes. This term does not find place 
in [s^ , since the calculations there have been carried out 
separately for the El and Ml transitions. If we identify 
2t:'^F^^ with 32A^/;s6 of [Hi where h denotes j, there 
is complete agreement between our expressions given by 
(|14|l and 115() for a and h and the corresponding expres- 
sions in j^j. If it is assumed that Els — Q and 

Sir" = Ell=^ = Ell-^ = Ely, (17) 

it follows that a, 6, c simplify to 

a = 8(|MUp + 3|Afl,n; 5 = 108|£;ii,p; c = 0, (18) 

leading to the beam analyzing power S(6') defined by (2) 
of 1331 which now assumes the form 



27 

E(0) - y|£;U|2sin2 0/i?, 



where the denominator 



27 

Z? = |Ml„p + 3|Ml,p + — |£;i„|2sin2| 



(19) 



(20) 



is proportional to the unpolarized differential cross sec- 
tion. The S(0) was determined experimentally a,t 9 = 
150° in ^ at 7-ray energy 3.58 MeV and at 61 = 90° in 
[S^l at seven 7— ray energies between 2.39 and 4.05 MeV. 
The measurements of S(6') in [s^ have led to empirical 
estimates of 



X = \M1VI\E1, 



(|MU|2-H3|Afl,n/|Sl,|2, (21) 



if Af Is is not set equal to zero. Under the same simplify- 
ing assumptions, it is interesting to note that the tensor 
target analyzing power is given by 



Al = ^[\Mly\^ - \\MU?]ID. 



(22) 



Thus experimental measurements of can lead to an 
empirical estimate of 



to study the energy dependence of R empirically in the 
energy region of interest to astrophysics. 

Finally, we may point out that, when the above sim- 
plifying assumptions are not made, the unpolarized dif- 
ferential cross section ^ itself is given by 



2^ 
6 



[a + bsm^ 9 — ccos^]. 



(25) 



where the coefficient c in third term can be determined 
by taking the difference between measurements of 
at two angles 9{^ '^1'^) and 71 — 9. It can also be 
determined in the same way from ^ given by (|13|l . For 
eg., Schreiber et al, 33] have measured ifT^ at 6 — 150° 
and for <j) — and 90°. Additional measurements at 
9 = 30° for the same angles (f> and at the same energy, 
could easily estimate c at 3.58 MeV. The measurements 
by Tornow et al., |34j at lower energies have been 
carried out at 6* = 90° and therefore not suitable for 
this purpose. It would therefore be desirable to carryout 
measurements at 6* 7^ 90° and at tt — at lower the 
energy, to determine c. The coefficient b is readily de- 
termined by taking the difference between (|13|l and H25|) 
at any angle 61 ^ or tt and for any value of </> ^ tt/A. 
Since b and c are thus known, one can determine a 
by measuring (|13|l or even H25|l . Thus a,b,c given by 
(I14|l . H15|) and Ijltill are determinable empirically without 
making simplifying assumptions as in |33l |. We may note 
from (|16|l that c goes to zero cither if Mlg is zero or if 
(|17|1 holds exactly. On the other hand, if an empirical 
determination leads to c 7^ 0, it implies simultaneously 
that Mis 7^ and the simplifying assumption (|17|l is 
invalid. 

Therefore an empirical determination of c appears de- 
sirable before carrying out the more incisive analysis of 
the experimental data suggested above. If c is found to 
be zero experimentally and H17|) is assumed to be valid, 
the measurements of (fHt along with ((T^ , (|^ and if^ 
hold promise for the more incisive empirical analysis, 
wherein R given by H24II also gets determined as a func- 
tion of energy along with (|21|l . where |Mlp represents 
\M l„p -I- 3|Mlsp. This will lead to a better understand- 
ing of the photodisintegration of deuterons at astrophys- 
ical energies of relevance for sharpening the predictions 
of BBN. 



Y^i\MU^--\Mls\')/\EU' 



(23) 



in the energy region of astrophysical interest. Since X 
and Y are known empirically as functions of energy, it is 
possible to estimate 



R = 



\Mls\^ ^ 2 (X-Y) 
|ML„|2 3{X + 2Y)'' 



(24) 
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